o 

(N 



(N 



Renormalized parameters and perturbation theory for an n-channel Anderson model 
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We extend the renormalized perturbation theory for the single impurity Anderson model to the 
n-channel model with a Hund's rule coupling, and show that the exact results for the spin, orbital 
and charge susceptibilities, as well as the leading low temperature dependence for the resistivity, are 
obtained by working to second order in the renormalized couplings. A universal relation is obtained 
between the renormalized parameters, independent of n, in the Kondo regime. An expression for the 
dynamic spin susceptibility is also derived by taking into account repeated quasiparticle scattering, 
which is asymptotically exact in the low frequency regime and satisfies the Korringa-Shiba relation. 
The renormalized parameters, including the renormalized Hund's rule coupling, are deduced from 
numerical renormalization group calculations for the model for the case n — 2. The results confirm 
explicitly the universal relations between the parameters in the Kondo regime. Using these results we 
evaluate the spin, orbital and charge susceptibilities, temperature dependence of the low temperature 
resistivity and dynamic spin susceptibility for the particle-hole symmetric n = 2 model. 
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I. INTRODUCTION 

The single impurity Anderson model^ has played an 
important role in understanding many aspects of the be- 
havior of electrons in systems with strong electron corre- 
lation. Non-perturbative methods have had to be devel- 
oped to make predictions for the behavior of the model in 
the strong interaction regime. Among the most success- 
ful have been the seminal and pioneering work of Wilson 
and associates 2 '? based on the numerical renormalization 
group (NRG), and the exact solutions using the Bethe 
ansatz for the linear dispersion version of the model^. 
Though the model was originally put forward to describe 
magnetic impurities in a host metal, it has proved to be 
applicable to other situations. One main area of applica- 
tion is as a model for strong correlation effects in quan- 
tum dots 6 . In this application certain parameters of the 
model, such as the impurity level which determines the 
electron occupancy on the quantum dot, can be varied by 
a gate voltage. This makes it possible to sweep through 
different parameter regimes of the model, which would 
be difficult to do for real magnetic inrpurities, and so the 
predictions of the model can be tested more rigorously. 
The presence of the narrow many-body resonance in the 
strong correlation (Kondo) regime at low temperatures 
can be inferred directly from the measurements of the 
current through the dot as a function of an applied bias 
voltage^. 

Apart from these direct applications of the model, it 
has also played a role in the calculations of strong correla- 
tion effects in lattice models. It is possible to map a class 
of infinite dimensional lattice models of strong electron 
correlation onto an effective Anderson impurity model 
with a self-consistency condition, which determines the 
density of states of the effective medium^.. This mapping 
requires that the self-energy is a function of frequency 
only which is the case in the limit of infinite dimension- 
ality, and the mapping is exact in this limit. For many 



strongly correlated systems it is known that the wave vec- 
tor dependence of the self-energy is much less important 
than the frequency dependence so this approach can be 
used as a good first approximation for systems in three 
dimensions (dynamical mean field theory (DMFT)). As 
the assumption of linear dispersion is not valid for the 
effective impurity model generated in this application, 
there are no exact Bethe ansatz solutions, so the most 
reliable non-perturbative approaches, such as the NRG 
have to be used. 

It has not proved possible so far to access the strong 
correlation regime of the Anderson model by an approach 
based purely on perturbation theory in powers of the lo- 
cal interaction U. However, it has been shown that, if 
the perturbation theory is reorganized such that the ba- 
sic parameters of the model are renormalized, then a per- 
turbation theory in the renormalized interaction U, taken 
only to second order gives formally the exact results for 
the low temperature properties and low frequency dy- 
namics, provided counter terms are taken into account 
to avoid over-countin g 10 ! 11 . The renormalized parame- 
ters have to be determined but these can be calculated 
very accurately from an analysis of the low energy exci- 
tations of an NRG calculation on the approach to the low 
energy fixed point 12 . So far this approach has only been 
developed in detail for the non-degenerate one channel 
model, but the approach is one that can be applied to 
a more general class of models including lattice models. 
Here we extend the calculations to an n-channel impu- 
rity Anderson model with the inclusion of a Hund's rule 
exchange term. The Hamiltonian takes the form, 
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where d^ , d m(T , are creation and annihilation opera- 



tors for an electron in an impurity state with total an- 
gular momentum quantum number I, and z-component 
to = —I,— I + 1, ...I, and spin component a =t>4-- 
The impurity level in a magnetic field H we take as 
tdma = £d - VbvH - fi B mH - /j,, where a = 1 (f) 
and c = — 1 (4-) and /i is the chemical potential, and 
/zb the Bohr magneton. The creation and annihilation 
operators c kma , c kma are for partial wave conduction 
electrons with energy e kma . The hybridization matrix 
element for impurity levels with the conduction electron 
states is V k ■ We denote the hybridization width factor by 
A mcr (e) = K^2, k \V k \ 2 S(e-e kmcr ), which we can take to be 
a constant A in the wide flat band limit. The remaining 
part of the Hamiltonian, lid describes the interaction be- 
tween the electrons in the impurity state, which we take 
to be of the form, 
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As well as the direct Coulomb interaction U between 
the electrons, we include a Hund's rule exchange term 
Jh between electrons in states with different to values. 
The sign for the exchange term has been chosen so that 
Jh > corresponds to a ferromagnetic interaction. This 
model can be used to describe transition metal impuri- 
ties, such as Mn or Fe, in a metallic host in the absence 
of spin orbit or crystal field splittings. We can interpret 
the model more generally with a = m + Z + Iasa channel 
index taking values a = 1, 2, ...n where n is the number of 
channels. The Hund's rule term tends to align the elec- 
trons on the impurity site such that for large U and large 
Jh the impurity state will correspond to a spin S = n/2. 
The model with Jh = has also been used to describe 
capacitively coupled double quantum dots^, where the 
impurity channels correspond to different dots. In that 
application, however, the inter-dot interaction U' will in 
general differ from the intra-dot interaction U, so the 
case here, with U' = U, is a special point with SU(2n) 
symmetry when Jh = 0. 

The structure of this paper will be as follows. In the 
next section we formulate the renormalized perturbation 
theory (RPT) for this model in terms of the renormal- 
ized parameters, id, A, U and Jh- We then show that 
the low temperature behavior, as measured by the charge 
and spin susceptibilities and the low temperature contri- 
bution to the resistivity, can be obtained exactly from 
the RPT taken to second order in powers of U and Jh- 
In the localized or Kondo regime we show that A, U and 
Jh can be expressed in terms of a single parameter which 
we take as the Kondo temperature Tk- This relation is 
independent of the channel index n and hence applies to 
all values of n. Though we cannot calculate A, U and Jh 
for the general n-model using the NRG we can calculate 
them for the two channel case n = 2. We look at this case 
in detail and confirm the universal relation between the 



renormalized parameters in the Kondo regime predicted 
using the RPT. 



II. RENORMALIZED PERTURBATION 
THEORY 



We start with the Fourier transform of the single par- 
ticle Green's function for the impurity d-state, 



Gd,a(oj) 
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{T T d ma (r)dl a (0))e^' T dr, (3) 



where uv = (2n/ + l)//3 and (3 — 1/T and the brackets 
(...) denote a thermal average. 
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where T, ma -(u) n , H) is the self-energy. For the zero tem- 
perature Green's function, which will be our main con- 
cern, uj n can be replaced by continuous variable u>, and 
summations over co n replaced by integrations over u>. For 
the perturbation theory in powers of U and Jh it will 
be convenient to separate the interaction terms in the 
Hamiltonian into the terms involving interactions be- 
tween electrons in the same channel and those between 
electrons in different channels. We rewrite the Hamilto- 
nian from Eq. ([2|) in the form, 

T~td = U 2, n d,m"[nd.mi 
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The vertices associated with the three types of interaction 
terms are illustrated in Fig. [TJ 
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FIG. 1: (Color online) The three interaction vertices corre- 
sponding to the terms in the Hamiltonian given in Eq. l|5|l. 



For the renormalized perturbation theory, the Green's 
function in Eq. (U) can be re-expressed as Gd t ma{^>n) — 



zGd t m(r(u n ), where G<2 imcr (w n ) is the quasiparticle 
Green's function given by 



Gd,ma\f^n) — 



1 



tuj n - l dma + ?Asgn(u/ n ) - T, m<T (oj n ,H) 

(6) 
and the renormalized parameters, tdma and A are given 
by 

e.dma = z(e d + £ mff (0, H)) - ^aH - [i B mH, A = zA, 

(7) 
where z = 1/(1 — 3E mo .(w, 0)/d(iui)) evaluated at ui = 0. 

The quasiparticle self-energy Y* mM (uj,H) is given by 



S md r(w,-ff) — 

2 S m(7 (u;, if) - £ mo .(0, ii) - uv 



diu> 



where we have assumed the Luttinger theorem 14 , 



ImS(O) = 0, so that ImS mCT (w) 



u> as w 



0. When 



expressed in this form, the uj = part of the self-energy 
and its derivative have been absorbed into renormaliz- 
ing the parameters tdma and A, so in setting up the 
perturbation expansion any further renormalization of 
these terms must be excluded, or it will result in over- 
counting. In working with the fully renormalized quasi- 
particles it is appropriate to use the renormalized or effec- 
tive interactions between the quasiparticles. In the single 
channel case, we defined the renormalized interaction U 
in terms of the four vertex T^^.f(ui,U2 1^3, W4) in the 
zero frequency limit 1 ^. In this case we need to consider 
the more general four vertex, r™^i™=£j(u>i,u;2, W3, W4), 
which corresponds to the Fourier coefficient of the con- 
nected skeleton diagram for the two particle Green's func- 
tion, 

{ T rd miai {T 1 )d m2a2 {T 2 )d\ n3a3 {T S )d mirTl {T i )), (8) 

with the external legs removed. Using the fact that 
the spin and angular momentum are conserved indepen- 
dently, and taking into account the antisymmetry condi- 
tions of the fermion creation and annihilation operators, 
it was shown by Yoshimorii^ that this vertex at zero fre- 
quency can be expressed in terms of two parameters, Tq 
and T P , as 



r mi<Ti;m 2 <T 2 fo 0) 
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To first order in the interaction terms, U and Jh, we have 
r<7 = U — Jh and T e — Jn. We generalize this result to 
specify the renormalized parameters, U , and Jh, by the 
relation. 
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where the factor z 2 arises from the rescaling of the fields 
to define the quasiparticle Green's function given in Eq. 
©. For n = 1 this reduces to 

z 1 r£;£(0, 0,0,0) = U{8%6% - 8%8%), (11) 

which is the definition of U used in earlier workiS. 

We can combine these terms to define a quasiparticle 
Hamiltonian if, 

" = / j e dmad ma d mcr + 2 J e kmaC kmrJ C kmtJ 
ma kma 

+ X! W ma c kma + V k *cl m J ma ) + H d (12) 
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(13) 
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The brackets :0: indicate that the operator O within 
the brackets must be normal ordered with respect to the 
ground state of the interacting system, which plays the 
role of the vacuum. This is because the interaction terms 
only come into play when more than one quasiparticle is 
created from the vacuum. 

The renormalized Hamiltonian is not equivalent to the 
original model, and the relation between the original and 
renormalized model is best expressed in the Lagrangian 
formulation, where frequency enters explicitly 11 . For 
simplicity, we consider the case in the absence of a mag- 
netic field, where the energy levels tdma are independent 
of to and a . If the Lagrangian density C(ed, A, U, Jh) de- 
scribes the original model, then by suitably re-arranging 
the terms we can write 

£{e dl A, U, Jh) = Cfa, A, U, J H ) + £ c (\i, A 2 , A3, A 4 ), 

(14) 
where the remainder part £ c (Ai, A2, A3, A4) is known as 
the counter term and takes the form, 

£ c (Al, A 2 ,A 3 ,A 4 ) = ^d m a{T)(X 2 d T - Xi)d m a(r) 
ma 

+ (A 3 -A 4 ) 22 d m ,a{T)dm'a'{T)dm'a'{T)dma(T) 

mm' era' 

+ A4 ^2 dma(T)d m 'a'(T)dm,a'(T)dm'a(T), (15) 

mm' aa' 

where Ai = -S(0), A 2 = z - 1, A3 = {z 2 U - U)/2 and 
A4 = (Jh — z 2 Jh)/2. Though we can express the coef- 
ficients Aj, i = 1,2,3,4, explicitly in terms of the self- 
energy terms and vertices at zero frequency, these rela- 
tions are not useful in carrying out the expansion. We 
want to work entirely with the renormalized parameters 
and carry out the expansion in powers of U and Jh- We 



assume that the A^ can be expressed in powers of U and 
Jh, and determine them order by order from the condi- 
tions that there should be no further renormalization of 
quantities taken to be already fully renormalized. These 
conditions are 



,(0,0) =0, 



dS mff (w,0) 



diuj 



0, 



(16) 



and that the renormalized 4-vertex at zero frequency, 
f^£™ (0,0, 0,0) is such that 



m3<73;m4(74 
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In the field theory context these conditions are more com- 
monly known as the renormalization conditions. They 
follow directly from the definitions of the renormalized 
self-energy in Eq. (|7|) and the definitions of the renor- 
malized parameters given in Eq. (J10[) . 

The propagator in the RPT is the free quasiparticle 
Green's function, 



G d,,na( UJ n) 



iuj n - e dma + iAsgn(uj n ) 



(18) 



The spectral density of the corresponding retarded 
Green's function gives the free quasiparticle density of 
states, p ma (uj) given by 



p^M = 



A/tt 



(W - e dm a) 2 + A 2 



(19) 



From Fermi liquid theory, the quasiparticle interaction 
terms do not contribute to the linear specific heat co- 
efficient 7 of the electrons. It follows that the impu- 
rity contribution to this coefficient is proportional to the 
free quasiparticle density of states evaluated at the Fermi 
level and is given by 



tE^(°)- 



(20) 



In the absence of a magnetic field this reduces to 7 = 
2n7r 2 p(°- ) (0)/3, where p(°\0) is the quasiparticle density 
of states per single spin and channel. 

If we integrate the free quasiparticle density of states 
in Eq. (JTSJ) to the Fermi level then we get (fidma) at 
T = 0, which is given by 



yldma/ 



TV 



1 



-tan 



(,!h 



A 



(21) 



which defines the phase shift rj ma in the channel with 
quantum numbers m and a. For this model it has 
been shown by Shiba 16 that (ridma) = Vma/^, giving a 
generalization of the Friedel sum rule, so that we have 



(ndma) = (ndrrur)] t ne quasiparticle occupation num- 
ber in each channel is equal to the impurity occupa- 
tion number in that channel. However, Yoshimori and 
Zawadowski 17 have shown that this form of the Friedel 
sum rule does not hold for a more general model in 
which scattering processes can occur between m-states, 
mi, 1712 — ► 7113,1714,, such that m\ + ra-i = mj, + m±. 
They derive a restricted form of the sum rule such 

that J2ma a mcr(ndmcr) = J2ma a maVd m a /^, where a ma = 

l,<j,m. In this more general case, therefore, the quasi- 
particle number does not equal the occupation number in 
the same channel but we have the more restricted result, 



E 



ma ""ma 



[ndma) = J2 m a a m<T(ndma)- Using either re- 
sult, however, we can derive expressions for the zero field 
spin Xs, orbital Xorb and charge \ c susceptibilities. We 
differentiate the combinations, ^2 m<y a- mc! (nd t ma), with 
Omur = Cj m an d 1 respectively, with respect to the mag- 
netic field or in the charge case with respect to id- To 
evaluate these expressions we need to calculate the renor- 
malized self-energy. This calculation taken to first order 
in U and Jh proceeds as in the one channel cas o 10 i 1:L , and 
gives 



Xs = 2n/4,5W(0)(l + {U + {n - l) Jh)p^(O)), (22) 



n(°)f 



Xorb 

and 



(n 2 -l)n%pM(0) 



12 



(l + (tf-3J H )p(°)(0)), 



(23) 



Xc = 2np^(0)(l-((2n-l)U-3(n-l)J H )pW(0)). (24) 

These results can also be obtained from a mean field 
theory on the quasiparticle part of the Hamiltonian given 
in Eq. (fl3|) 18 ' 19 . It can be shown using the Ward identi- 
ties derived by Yoshimori 15 , which are generalizations of 
the Ward identities derived by Yamad a 20 ' 21 for the sin- 
gle channel case, that these results are exact. Hence all 
higher order correction terms in U and Jh cancel out. 

In the localized regime a large value of U suppresses the 
charge fluctuations on the impurity so Xc ~ 0. Treating 
this as an equality, we get a relation between p(°'(0), U 
and Jh, 



((2n-l)U- 3(n - 1) J H )p (0) (0) = 1 
When Jh = 0, this reduces to 

Up (o) {0) 



(2n- 1)' 



(25) 



(26) 



For the case of half-filling, where id = and p^ (0) = 
1/ttA, the non-linear relation between the renormalized 
parameters in Eq. ([23)) becomes a linear relation between 
A, U and Jh, 



ttA = (2n - 1)U - 3(n - 1) J H - 



(27) 



For Jh = we get 



U 



irA 



(2n-l)' 



(28) 



which agrees with the one channel result U — irA for 
n = 1. 

When Jh = and we are in the localized limit, we 
have only one energy scale which we can take to be 
the Kondo temperature, defined for general n such that 
7 = 7r 2 n/6T K , equivalent to taking /5 (0) (0) = 1/4T K . 
In this limit the result for the Wilson ratio, Ryf = 
7r 2 x s /3/i B 7 = 2n/(2n — 1). This is the same as that for 
the N-fold degenerate Anderson model used to describe 
rare earth impurities for N = In. This result could have 
been anticipated, because the models can be shown to 
be equivalent by putting the orbital m and spin indices 
a into a combined index v = (m, c) 2 ^. 

Switching on the interaction Jh (> 0) will reduce the 
local orbital fluctuations, as the configuration with the 
spins aligned will be favored. For Jh 3> ttA we can ex- 
pect the orbital fluctuations to be almost fully suppressed 
so that Xorb ~ which, as an equality, gives a further re- 
lation between p(°)(0), U and Jh, 



(3J h -C/)p (0) (0) = 1. 



(29) 



At half-filling this gives another linear relation between 
A, U and Jh, 



ttA = 3 J H - U. 



(30) 



An equivalent condition to that in Eq. (|30[) can be ob- 
tained using the argument of Nozieres and Blandin 23 that 
the occupation number in a channel m should be inde- 
pendent of any small change in the chemical potential 
in a channel m! ^ m in this regime. When both the 
local charge and orbital fluctuations are suppressed, the 
renormalized parameters can be expressed in terms of 
the Kondo temperature Ik- From Eq. (|25|) and ([29]) we 
deduce 



U = | J H = 47k. 



(31) 



for the particle-hole symmetric case we have l/p(°\0) = 
7rA = 4Xk, so then we have 



ttA = U =-Jn = 4T K . 



(32) 



which was conjectured earlier on the basis of a phe- 
nomenological mean field approac h 18 ' 19 . A notable fea- 
ture of this result is that there is no explicit dependence 
on n. In this regime from Eq. (|22p we have for the spin 
susceptibility, 



where S — njl and g - 
R w = 2(n + 2)/3i^3. 



(gu B ) 2 S(S+l) 



3T K 



(33) 



2, which leads to a Wilson ratio, 



Yoshimorii^ has also derived an exact result for the 
low temperature impurity contribution to the resistivity 
in the particle- hole symmetric case and H = 0. In terms 
of the renormalized parameters, the result is 



R{T) =Ro[l 



ir i (l + 2I R )T 2 



0(T 4 ) 



(34) 



where Ir is given by 



lR = (~p {0 \0)) 2 {{2n-l)U 2 -Q{n-l)MU-Jn))- (35) 

This result can be derived in the RPT from a calculation 
of the renormalized self-energy S(w) to second order in U 
and Jh- With the Hund's rule interaction term, there are 
several types of second order scattering diagrams which 
are illustrated in Fig. [5] The vertices are of the same type 
as shown in Fig. Q] but are weighted by the renormalized 
interaction terms. The calculations follow along similar 
lines to those for the single channel case n = l 10 ' 11 . The 
first order diagrams and the terms linear in uj are can- 
celed by the counter terms to this order, and there are no 
corrections from the counter terms to the vertices to sec- 
ond order for the case with particle-hole symmetry. The 
contributions to Ir from diagrams of the types (i) to (iv) 
respectively in units of (/^^(O)) 2 are: U 2 ; 2(n — 1)J H ; 
2(n - 1)(U - Jh) 2 ; -2(n - 1)J H (17 - Jh); which give the 
result in Eq. (|33| . 






(iv) 



FIG. 2: (Color online) Second order diagrams in the renor- 
malized perturbation theory. 



In the localized regime at half-filling the result in Eq. 
simplifies to give 



R{T) = Ro[l 



7T 4 (5 + 4n) ( T\ 2 



')(> 



Tk 



+ 0(T 4 ) , (36) 



which agrees with the result derived by Nozieres^ 4 . and 
Yamada 20 ' 2 ^ for the case n = 1. Thus all the exact Fermi 
liquid relations can be derived from the RPT taken to 
second order only. 

It was shown in earlier work 2 ^ that the RPT approach 
can provide a description of the dynamic spin susceptibil- 
ity for the n = 1 model in the low frequency regime. The 



calculation takes account of the repeated quasiparticle 
scattering, giving results which are exact in the low fre- 
quency limit uj — > 0, and in remarkably good agreement 
with the results from a direct NRG calculation. We ex- 
tend the calculation to the n-channel model given in Eq. 
([I]) and ([2]). We consider the Fourier transform of the 
transverse spin susceptibility, 



which leads to the result, 



X. 



+- 



m m' 

(37) 



where 



= 2im'/l3 and S 



d,m 



^mi^mt \^d, 



— rfi rl 



s. 



d,m 



{ n d,-\ — "d,4.)/2). We consider the scat- 
tering of a spin up quasiparticle with a spin down quasi- 
hole both in channel to, in the absence of a magnetic 
field. This particle-hole pair can scatter into a particle- 
hole pair in the same channel to or a different channel 
to' 7^ m. We consider the scattering into the same chan- 
nel first of all. The matrix element for this process is U, 
except we must allow for the fact that U already takes 
into account these processes for uj — so, to prevent 
over-counting, we must use U — A3, where A3 is the cor- 
responding counter term. It will be convenient to use the 
notation U for U — A3. Just taking this type of repeated 
scattering into account gives us a result which has the 
same form as in the single channel case n = 1—, 



xt ( w + iS ) = 4^/4 



IL+-(u + i6) 

l-UU+-(uj + iS) : 



(38) 



where we have analytically continued to real frequency 
ui. The free quasiparticle-quasihole propagator in a single 
channel, II -1 (u+iS), is independent of the channel index 
in the absence of a magnetic field, and is given by 



U + -{u + i8) = 



A 



ttuj(uj + 2iA) 



n(? d + A*y 



In 1 



w = 



£ d + iA 



In 1 - 



e d -»A 



w^O, (39) 



for S — > +0. We must also take into account that the 
quasiparticle-quasihole pair being created in channel to 
can scatter into a different channel to', and also be fi- 
nally annihilated in a channel with to' ^ to. The ma- 
trix element for this type of scattering is Jh , correspond- 
ing to the diagram in Fig. Q] (ii), but again, to avoid 
over-counting, we replace it by Jh- In the absence of a 
magnetic field, the quasiparticle-quasihole propagator is 
independent of the channel index to, so the summation 
over the states to' introduces a factor n — 1. The result 
of taking these scattering processes into account is that 
the pair propagator II H (w + iS) in Eq. (|38|) is replaced 
by 



IL+-(uj + i5) 



J H (n-l)n+-(uj + i5) 



(40) 



Xt (0J + i8) 



4timb: 



IL+-{cu + i5) 



' 1 - (U + (n - 1) Jff )n+- (w + iS) 

(41) 

We need to determine the combination U + (n — 1) Jh- 
We can do this by requiring that this expression gives 
2\ s in the zero frequency limit, which is equivalent to the 
requirement that these scattering processes contribute to 
the four vertex at zero frequency are not over-counted. 
This condition gives 



U + {n- 1) J H 



U + (n-l)J H 



l + ([7 + (n-l)J H )p(°)(0) 



(42) 



In the Kondo regime this condition simplifies to U + (n — 
\)Jh = 2TK(l + 2n)/(2 + n), which gives the one channel 
result U = 2T K for n = 1. 

By rewriting Eq. (|41[) in the form, 



4n^L 



X t~(u+iS) IL+-(uj + iS) 



-(U+(n-l)J H ), (43) 



and taking the imaginary part, it is straight forward to 
show that the expression for x s (w) satisfies the exact 
Korringa-Shiba relation, 



lim 

u-10 



Imx + (w + id) _ nx s 



nul' 



(44) 



which was proved for this model by Shiba 16 and more 
generally by Yoshimori and Zawadowsiii. 

So far we have not discussed how one can calculate 
the renormalized parameters e^, A, U and Jh- In the 
Kondo regime these reduce to a single parameter Tk, so 
one possibility is to deduce its value from experiment by 
fitting the predictions to the measurements of a physical 
quantity in the low temperature regime, say the impurity 
susceptibility or resistivity. Outside the Kondo regime we 
have four parameters to determine, and to calculate all 
four from experiment one loses much of the predictive 
power of the approach. However, it was shown earlier for 
the single channel Anderson model how the parameters, 
Cd, A and U, can be calculated in terms of the bare pa- 
rameters, ed, A and U, from the many-body low energy 
excitations of an NRG calculation 12 . There are problems 
in carrying out this procedure for the general n-channel 
model, due to the truncation of states which has to be 
carried out in an NRG calculation to reach the very low 
energy scales. Truncation means that only a fraction 
1/4™ states can be retained at each NRG iteration. It is 
possible, however, for the case n = 2 to compensate for 
the lower percentage by increasing the number of states 
kept at each iteration as the matrices do not get so large. 
In the next section we present for calculations of A, U, 
and Jh, in terms of A, U, and Jh, for the n = 2 model. 



III. NRG CALCULATION OF THE 
RENORMALIZED PARAMETERS FOR N=2 

The two-channel model the Hamiltonian T-Ld given in 
Eq. ([2]) can be re-expressed in the form, 



u d = u Y, 



n<iat n dal- 



■Ui 2 '^2n di i a n di 2a>- 



-2JnSd,i-Sd,2, 



(45) 
with a ferromagnetic Heisenberg exchange coupling 2 Jjj 
between the electrons in the different channels, and U12 — 
U — 3Jh/2. Our calculations will be restricted to the 
particle-hole symmetric model so we take td — —U/2 — 
U\2 in the one-electron part of the Hamiltonian given in 
Eq. |1|. The energy of the two electron triplet state 
of the isolated impurity with particle-hole symmetry is 
— 2U + Jh and that of the 4-electron or 0-electron state is 
0, so if we are interested in the case when the triplet state 
is the ground state configuration, we need to consider the 
regime U > Jh/2. 

For the NRG calculations the model is recast in a form 
such that the impurity is coupled via a hybridization V to 
two tight binding chains which describe the conduction 
electron states, one chain for each channel. The con- 
duction electron band is discretized with a discretization 
parameter A > 1, such that the couplings decrease along 
the chains as A~ N/>2 for large N, where N is the iVth 
site along the chain from the impurity. The calculations 
are then carried out iteratively by direct diagonalization, 
starting at the impurity site and adding one further site 
to each chain at each iteration step. The number of basis 
states used has to be truncated when the matrices get too 
large for diagonalization on a practical timescale, which 
can occur after only a few iterative steps. When trun- 
cation is applied a fixed number of states is retained at 
each step. For the n — 2 model considered here, we take 
3600 states, which is a factor of 3 to 4 more than for the 
non-degenerate model (n = 1) and a discretization factor 
A = 6. We can check the expected accuracy of our calcu- 
lations by using this value for A to calculate U and wA 
for the single channel model and compare with the values 
deduced indirectly from the exact Bethe ansatz results 
for the specific heat coefficient 7 and the zero tempera- 
ture spin susceptibility-^. For U/nA — 2, 7rA = 0.01, 
keeping 900 states, we get the values, U — 0.2295 and 
7rA = 0.2387, which can be compared with those deduced 
from the Bethe ansatz, U = 0.2301 and ttA = 0.2392. 
This gives an accuracy of better than 0.3%. For further 
details on setting up the NRG calculations, we refer to 
the original papers^ and the recent review article^. 

With this discrete spectrum the Green's function in 
Eq. ((4]) takes the form, 



GdA") 



1 



iu - e dma - \V\ 2 g aa (iu) - S mCT (w) : 



(46) 



The connection between the NRG approach and the 
renormalized perturbation theory is based on identify- 
ing the quasiparticle Hamiltonian, given in Eq. (fT2|) and 
(fT"3")l , as the low energy fixed point of the NRG together 
with the leading irrelevant terms^. The lowest single- 
particle excitations from the NRG ground state should 
correspond to a quasiparticle excitation described by the 
one-body part of the quasiparticle Hamiltonian as given 
in Eq. (|T2"]) . For the calculation of the interaction terms, 
U and Jh , from the NRG we have to consider the differ- 
ence between two-body excitations from the NRG ground 
state and the two corresponding one-body excitations. 

The low energy single-particle excitations are given 
by the poles of the non-interacting quasiparticle Green's 
function when analytically continued to real frequency uj. 
The equation for these poles is the same as that for the 
non-interacting model but with a renormalized hybridiza- 
tion V and energy level id- Therefore, the lowest energy 
single particle and hole excitations, E p (N) and Eh(N), 
from the interacting ground state should be solutions of 
the equation, 



uj-e d -\V\ 2 9a(7 (Lo)=0. 



(47) 



where g atT (iuj) is the Green's function for the first site for 
the isolated conduction band chain. 



If we substitute the excitations energies, E p (N) and 
Eh(N), as calculated in the NRG for a finite chain length 
N, into Eq. (|47| then we can deduce corresponding 
Af-dependent renormalized parameters U(A^) and £d(N). 
Only if V(N) and £d{N) become independent of N for 
large N, do the low energy one-particle energy levels of 
the interacting system correspond to those of a renor- 
malized non-interacting model. If this is the case, then 
the asymptotic values for large A~ define the renormalized 
parameters V (and hence A) and e d - 

To calculate the renormalized interaction terms, we 
first have to diagonalize the non-interacting impurity 
model with the renormalized parameters, which describes 
the quasiparticles. The interaction terms are then added 
to the quasiparticle Hamiltonian and expressed using the 
diagonalized single quasiparticle states as a basis. The 
energy difference between the lowest two-particle state 
and the sum of the corresponding two quasiparticle states 
is equal to the expectation value of interaction terms 
in the quasiparticle Hamiltonian. The interaction pa- 
rameter U can be calculated from the NRG results for 
the lowest two-particle excitation in the same channel 
which will be independent of Jh- For a finite length 
chain N the value U(N) will depend upon N, and for 
this to correspond to a low energy quasiparticle Hamil- 
tonian U(N) should become independent of N for large 
N. The asymptotic values of U(N) for large A~ defines 
the renormalized parameter U. Similarly, to calculate Jh 
we look at the difference between the single and triplet 
states of a two-particle excitation with one electron ex- 
citation in each of the two channels. This excitation will 
be independent of U and depend only on Jh- Using the 
NRG results for a finite chain of A" sites, we can define a 
parameter Jn(A r ), with Jh given by the asymptotic value 



of Jn(N) for large N. Further details on the calculations 
of the renormalized parameters from the low energy NRG 
states can be found in reference 12 . 

We first show results for the renormalized parameters 
as a function of N. We show a typical case in Fig. [3] for 
the parameters U(N), ttA(N) and 3Jh(^V)/2 as a func- 
tion of N for vrA = 0.01, U/nA = 3.6 and J h /ttA = 0.15, 
which is a parameter set corresponding to a point in the 
Kondo regime where the orbital fluctuations have been 
suppressed. The results demonstrate that not only is 
there a plateau region for all the parameters for large 
N, but also that the asymptotic values of U(N), ttA(N) 
and 3Jh(-/V)/2 correspond to a single energy scale and 
satisfy the relation given in Eq. ([32| . The choice of a 
relatively large value of A = 6 means that the conver- 
gence to a plateau region is achieved for relatively small 
values of N. The plateau region is finite because the 
renormalized parameters correspond to the leading irrel- 
evant corrections to the free fermion fixed point of the 
Wilson renormalization group transformation 2 , so even- 
tually they diverge from the plateau when N is such that 
the decreasing irrelevant corrections become of the same 
order as the uncertainties in the numerical computation. 
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FIG. 3: (Color online) A plot of A(JV)/ttA, U(N)/tyA and 
3J H (iV)/27rA versus N for U/ttA = 3.6 Jh/ttA = 0.15 and 
7rA = 0.01. The inset shows the convergence of these param- 
eter to a common limit in this case as the bare parameters 
correspond to a point in the Kondo regime. 



We next look at the renormalized parameters in the 
different parameter regimes of the model. In Fig. |4] we 
show the results for A/3A and U/ttA versus U/irA for 
J H = (ttA = 0.01). We predicted from Eq. ^ that 
for large U/ttA we should have a single energy scale such 
that for n = 2, U = ttA/3 and the results clearly show 
that this is the case for U/irA > 3. The numerical results 
for the ratio U/ttA for large U give the value 1/3 to an 



accuracy of 0.01%. 

In Fig. [5] and [5] we compare the results for these two 
quantities with those for the single channel model n = 1. 
We can see that the parameters A and U the fall off 
with increase of U much more slowly for the two chan- 
nel model. This is because in the two channel model we 
have unsuppressed fluctuations of the orbital component. 
When Jh = and finite U, at half-filling in the isolated 
impurity for the two channel model there are six degen- 
erate two-electron configurations with energy 2e c [ + U. 
Both the n = 1 and n = 2 models in the Kondo regime 
can be described by localized SU(2n) Kondo model. For 
the case n = 1 it is the s-d or SU(2) Kondo model and 
for n — 2 the Coqblin-Schrieffer or SU(4) Kondo model. 
The Hamiltonian for the SU(2n) Kondo model takes the 
form, 

H K (2n) = Jeff Y^ Y »y c t'y c k, v + Y e k4,* ( V> 



v.v'k.k' 



i/.k 



(48) 
where the sum over v = 1, 2, ...2n, and with particle -hole 
symmetry J c s = A\V\ 2 /U. The operators Y v y obey the 
SU(2n) commutation relations, 



\Y , Y ,, ;»1 



Y II iS n 



(49) 



with ^2i V Y v>v = nl. For n = 1, Y v y = \v){v'\, where 
\v) are the single electron impurity states with spin up 
(y = 1) and spin down (y — 2), giving a two dimensional 
representation for the Y u y . In the two channel case for 
half-filling the representation of the operators Y v y is six 
dimensional and details of the Y v y in terms of the two 
electron impurity states are given in the Appendix. 

The relation Tk = 7rA/4 applies for both the n = 1 
and n = 2 models in the case of particle-hole symme- 
try. In the single channel case n = 1, Tk is known from 
the Bethe ansatz solution, and is given by Tk/ttA — 
v /^/2^ e -T 2 «/8+0.5/« ) where u = u/ n A 2S 7 and the NRG 

results for Tk deduced from 7rA are in precise agreement 
with this expression for large U. For n > 2 there is 
no Bethe ansatz solution for the model with finite U. 
However, there is a Bethe ansatz solution for the SU(N) 
Kondo model (Coqblin-Schrieffer model) and the N-fold 
degenerate Anderson model with U = od^r— which gives 
in the exponential for Tk a factor proportional to 1/N . 
The prefactor is not universal and depends on the cut- 
offs used for the high energy excitations in the model. 
We have taken for the two channel case, therefore, the 
expression Tk/tt A = 1.0lMe- 7r2 "/ 16+0 - 25 / u /27r, where the 
prefactor has been chosen to give the most reasonable fit 
to the data. The result of this fitting is shown in Fig. [71 
where it can been seen that the agreement is very good in 
the strong coupling range U/ttA > 4.0. The same form 
for Tk was used in reference^ 2 ., and found to be in good 
agreement with their NRG results. 
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FIG. 6: (Color online) A comparison of U/nA versus U/nA 
for the n — 1 and n = 2 models for Jh = and it A — 0.01 



FIG. 4: (Color online) A plot of A/37rA and U/nA versus 
U/nA for J H = and ttA = 0.01. 



this regime the relations between the renormalized pa- 
rameters are such that U12 =U — 3Jh/2 = 0. 
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FIG. 5: (Color online) A comparison of A/A versus U/nA 
for the n = 1 and n = 2 models for Jh = and 7rA = 0.01. 
The inset showsthe corresponding values for the Wilson ratio 
Rw = 1 + U/nA. 
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FIG. 7: (Color online) A plot of T K /nA (ttA = 4T K ) as a func- 
tion of U/nA for Jh = and 7rA = 0.01. The dashed curve 

corresponds to the formula 1.01ue~' r2u/16+0 - 25/ " /2tt, where 
it = U/nA. The inset shows a plot of the logarithm for the 
same two curves. 



In Fig. [8]we look at the effect of switching on the Hund's 
rule term Jh for a relatively large value of U, U/nA = 
4.0, which is sufficient to suppress the charge fluctuations. 
As we increase Jh, we begin to suppress also the orbital 
fluctuations, such that when Jh/7tA > 0.1 we are in the 
regime where we have a single energy scale. This we refer 
to as the Kondo regime with the Kondo temperature in 
the particle- hole symmetric case given by 7rA = 4Tk. In 



In Fig. [5] we plot the corresponding spin, orbital and 
charge susceptibilities using the expression for these given 
in Eq. I|22|l , (|23|) and l|24|l for the set of parameters used 
for Fig. [5) The fact that the charge susceptibility is 
almost zero, due to the large value of U, U/nA = 4.0, 
means that the renormalized parameters must satisfy Eq. 
([2"T]) . This provides some insight into why the value of U 
increases initially as Jh is switched on. For small Jh the 
change in Jh is almost linear whereas the change in 7rA 
is relatively small. Therefore to satisfy Eq. (|2"T)) U must 
also increase almost linearly in this region. We can also 
see from Fig. [5] that the orbital susceptibility is small 



1.0 



(multiplied by a factor 10 in the figure), and decreases 
monotonically as Jh increases. 





FIG. 9: (Color online) A plot of the spin susceptibility Xs 
(units of 4/ig), 40 x Xc where Xc is the charge susceptibility, 
and the 10 xxorb, where Xorb is the orbital susceptibility (units 
of jUg/4), versus Jh/7tA for the same parameter set as in Fig. 

El 



FIG. 8: (Color online) A plot of ttA/ttA, U/nA 3Jh/2ttA 
versus Jh/wA for U/ttA = 4.0 and A = 0.01. There is a single 
renormalized energy scale when Jh/ttA > 0.1. The inset 
shows the corresponding Wilson ratio Rw — 1 + (U+ Jh)/ttA. 



In Fig. [TU]we explore a different parameter regime. Here 
the parameters A, U, and 3Jh/2 are plotted for a range 
of values of U for Jh/VA ~ 0.05. We see that for this 
smaller value of Jh a large value of U/ttA ~ 5.5 is re- 
quired before the orbital fluctuations are suppressed and 
the Kondo regime is achieved. We suggest that the ex- 
planation for this behavior is that a large U strongly 
renormalizes the effective hopping parameter, which is 

proportional to V A, so that the relatively weaker Jh is 
then sufficient to suppress the orbital fluctuations. The 
inset of Fig. [TU] shows that the renormalized parame- 
ters do actually converge for U/ttA > 5.5 to a com- 
mon value. In Fig. [11] we plot the Wilson Xs/l ratio, 
Rw = 1 + {U + Jh)/7tA, for the parameter set given in 
Fig. [TO] It shows a steady increase from a value i?vv ~ 1 
for small U with a leveling off at U/irA ~ 5 and then 
a convergence to the value i?w = 8/3, corresponding to 
that of the localized S = 1 two channel Kondo model. 



In the Kondo regime for the model with Jh ^ with 
particle-hole symmetry we have Tk = vtA/4. This regime 
occurs when Jh is large enough so that the triplet state 
of the impurity has a much lower energy than the other 
two-particle impurity states. The effective coupling of 
this state to the conduction electrons, via virtual transi- 
tions to either single particle or three particle impurity 
states induced by the hybridization, leads to an exchange 



h \ 



0.8 



0.6 



0.4 



- V 



0.2 



* 7lA/7lA 

X 



N, 



s* 



****** 




4 4.5 5 5.5 6 6.5 7 

"*^ U/nA 

U/nA *%^ 
<***************^. *% % 



'^ri^*********^^ 




12 3 4 

U/tuA 

FIG. 10: (Color online) A plot of ttA/ttA, U/tvA and 
3J H /27rA versus U/tyA for Jh/ttA = 0.05 and A = 0.01. 
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FIG. 11: (Color online) A plot of the Wilson ratio Rw — 1 + 
(U + Jn)/nA versus U/ttA for Jh/ttA = 0.05 and A = 0.01. 
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model of a localized spin 1 coupled to the two channels 
of conduction electrons with an effective antiferromag- 
netic exchange interaction J c g = AV 2 /(U + Jh)- This 
in turn will lead to a Ju-dependent term in the Kondo 
temperature of the form Tk ~ exp(— an 2 Jh/7tA), where 
a is a dimensionless numerical coefficient. This implies 
that in the Kondo regime Tk will vary exponentially with 
Jh/ttA. In Fig. [TJ we plot T K from the NRG results 
against Jh/ttA and compare them with an exponential 
fit. The inset shows the plot of the logarithm of Tk, 
Ln(TK/7rA), versus Jh/ttA. It can be seen that the expo- 
nential form does fit well with the results for the Kondo 
range Jh/ttA > 0.1 with the value a = 1.49. There 
is a slight deviation for the largest values of Jh shown, 
but the coefficient a depends on the range chosen for the 
curve fitting. 
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FIG. 13: (Color online) A plot of the real part of the dynamic 
spin susceptibility \s (w) (in units of 8/i B ) f° r Jh/ttA = 0.05 
and Jh/ttA = 0.15 with U/ttA = 4, A = 0.01. 
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FIG. 12: (Color online) A plot Tk/ttA (full curve) and 
0.0854exp(-1.49Tr 2 Jh/ttA) (dotted curve) versus Jh/ttA for 
U/ttA - 4, A = 0.01. The inset shows the Ln(T K /TrA) and 
-1.49tt 2 Jh/ttA + Ln(0.0854). 
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FIG. 14: (Color online) A plot of the imaginary part of 
the dynamic spin susceptibility Xa~ _ ( w ) ( m units of 8/i B ) for 
Jh/ttA = 0.05 and Jh/ttA = 0.15 with U/tvA = 4, A = 0.01. 



IV. CONCLUSIONS 



Using the renormalized parameters for values of 
Jh/ttA = 0.05 and J h /ttA = 0.15 taken from fig. M 
corresponding to U/irA = 4 and A = 0.01, we have eval- 
uated the expressions for the dynamic spin susceptibility 
given in Eq. (|4"Tj) . The result for the real part is shown 
in Fig. Q2J It illustrates the narrowing and height in- 
crease of the central peak with the larger value of Jh- 
In Fig. [TJthe imaginary part of xt ( w ) is shown. The 
marked increase in the change of the gradient through 
the origin for the larger value of Jh, can be explained 
as a consequence of the Korringa-Shiba relation given in 
Eq. (0). 



The point of this study here for the n-channel Ander- 
son model has been to show how the renormalized per- 
turbation approach (RPT) can provide an asymptotically 
exact way of calculating the low temperature and low fre- 
quency behavior of the model in all parameter regimes. 
There have been many previous studies of related multi- 
orbital impurity models using a variety of approaches. 
The general n-channel Anderson model with finite U 
has not so far been solved using the Bethe ansatz, but 
there are exact solutions using this technique for the n- 
channel Kondo model coupled to a spin 6 * 29 i 30 ' 32 and the 
n-channel Anderson model in the infinite U limit. In the 
latter case the impurity occupation number is restricted 
to the range rid < 1— • The main focus of the work on 
the n-channel Kondo model, however, has been on the 
over-screened case for n > 25, where 5 is the spin of 
the impurity, as the model has a low energy non-Fermi 
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liquid fixed point. There have been many NRG studies 
of multi-orbital models and this work has been surveyed 
in the NRG review article^. The concern in most of the 
NRG work, however, has been with the calculation of the 
one-electron spectral densities, and mainly for the mod- 
els without the Hund's rule term. There has also been 
a recent study for the Jh = model using the local mo- 
ment approach, which includes NRG calculations for the 
case n = 222, and NRG studies of capacitively coupled 
quantum dot*^. 

The main feature of the RPT approach is that the 
calculations are carried out in terms of renormalized pa- 
rameters which have a clear physical meaning in terms 
of the quasiparticles and their interactions. For the n- 
channel model they correspond to renormalizations of the 
parameters, e ( i, A, U and Jh, which specify the model. 
In the strong correlation or Kondo regime all these pa- 
rameters can be determined explicitly in terms of a single 
low energy parameter, the Kondo temperature Tk- For 
the case n = 2, we have been able to deduce the renor- 
malized parameters from the low lying excitations in an 
NRG calculation. The NRG results have confirmed the 
relations we derived between the renormalized parame- 
ters in the Kondo regime. As we have explicit expressions 
in Eq. ([50" ]) . (|2"2" ]l . ([25] ) and (J2J|| for the specific heat co- 
efficient, spin, orbital and charge susceptibilities at zero 
temperature, these quantities were calculated simply by 
substituting the renormalized parameters into the rele- 
vant formulae. This procedure is very accurate and by- 
passes the usual NRG method which involves a subtrac- 
tion procedure to isolate the impurity component. As 
there is a large parameter space to explore we have re- 
stricted the NRG calculations here to the particle-hole 
symmetric case. However, the RPT results are valid in 
all parameter regimes and the behavior of the model away 
from particle-hole symmetry will be the subject of a sep- 
arate publication. 

In setting up the RPT no approximation has been 
made, other than the assumption that the self-energy 
and its derivative are real and non-divergent at the Fermi 
level u — 0. This means that there is the possibility 
of extending the results to higher temperatures and fre- 
quencies. Some preliminary results have been achieved 
by including diagrams beyond second orde r 11 ' 33 for the 
single channel model and this topic is currently being 
studied. The RPT in the Keldysh formalism can also be 
applied to non-equilibrium behavior and has been applied 
to the calculation of the non-linear corrections to the dif- 
ferential conductance for a quantum do t 34 i 35 , including 
an arbitrary magnetic field 3 ^. 

The RPT approach is not restricted to impurity mod- 
els, but the calculation of the renormalized parameters 



for lattice models presents more of a problem as the 
NRG method cannot in general be applied. However, 
for infinite dimensional lattice models one can use the 
dynamical mean field theory (DMFT) to map the model 
into an effective impurity one, so the NRG method can 
then be used. This approach has been used to calculate 
renormalized parameters for the one-band Hubbard and 
Hubbard-Holstein model o 37 ! 38 . The work presented here 
opens up the possibility of extending this method to the 
two-band Hubbard model with a Hund's rule coupling. 
We have found that the Hund's rule term plays an impor- 
tant role in enhancing the magnetic response in the 2-fold 
degenerate model. It is known that the single band Hub- 
bard model does not provide a basis for explaining the 
occurrence of ferromagnetism in 3d metals, as it predicts 
a ferromagnetic ground state only in a very restricted pa- 
rameter regime, very close to half-filling and for a value 
of U much greater than the band width. It is likely that 
the inclusion of the Hund's rule coupling is essential to 
describe ferromagnetism in 3d materials. 
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Appendix 

For the particle- hole symmetric model with n = 2, the 
model given in Eq. (|48p can be derived by taking account 
to order \V\ 2 the effects of virtual excitations from the 
2-electron to the local 1-electron and 3-electron impurity 
states. We denote the 1-electron basis states, |1 f), |1 4), 
|2 t) an d |2 |), by \v) with v = 1,2,3,4 respectively. 
The 2-electron states we denote by \u,v'), with v ^ v' 
and \v\v) represents the same state. This gives a six 
dimensional basis set. In terms of the Hubbard operators 
X(y,v>):{v», v »>) = \v,i/){v",v"'\, the Y„y are given by 



Y„ 



J2 (-ir*(^»):(">")> (50) 



v" ^v,v" ^v' 



for v' > i/, where a = 1 if v < v tl < v' . otherwise a = 0. 
The Y v i ^ for v f > v can be deduced from (|5Q)) using 
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